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Abstract 

We investigate the emission of gravitons by a cosmological brane 
into an Anti de Sitter five-dimensional bulk spacetime. We focus on 
the distribution of gravitons in the bulk and the associated produc- 
tion of 'dark radiation' in this process. In order to evaluate precisely 
the amount of dark radiation in the late low-energy regime, corre- 
sponding to standard cosmology, we study numerically the emission, 
propagation and bouncing off the brane of bulk gravitons. 
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1 Introduction 

In the last few years, a new cosmological scenario has emerged, based on 
the assumption that our universe is a brane: a sub-space embedded in a 
larger, bulk space, with additional dimensions. In this context, a model 
that has attracted particular attention is that of a self-gravitating three- 
brane embedded in an empty five- dimensional spacetime PP, and especially its 
extension j2j 0] inspired by the (non cosmological) Randall-Sundrum model 

® 
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In the latter case, where the bulk is endowed with a negative cosmolog- 
ical constant, it is possible to find a viable cosmological model, compatible 
so far with the available data. Due to the cosmological symmetries, it can 
be shown that in this model the most general bulk geometry corresponds to 
a portion of five-dimensional Anti de Sitter-Schwarzschild (AdS-Sch) space- 
time, described by the metric 

ds 2 = -(k + //V - ^) dt 2 +(k + /iV - 1 dr 2 + r 2 dY? k , (1) 

with k — — 1, 0, 1 depending on the curvature of the three-dimensional spa- 
tial slices. The bulk cosmological constant is related to the mass scale \i 
via 

A = -6/i 2 . (2) 

From the point of view of the brane, whose energy density is supposed 
to be the sum of an intrinsic tension a and of the usual cosmological matter 
energy density p, cosmology is governed by the brane Friedmann equation, 
which is different from the standard one and reads [3] 

This result can be obtained by applying the gravitational junction conditions 
for a moving brane in (Jlj), the cosmological scale factor a(t) being simply the 
radial coordinate r at the brane position. The last term in eq. (|5|L C/a 4 , 
is usually called "dark radiation" or "Weyl radiation" and represents the 
influence of the bulk geometry on the cosmology in the brane. 

If the bulk is strictly AdS-Sch then C is a constant and corresponds to 
the five-dimensional mass in the Schwarzschild AdS metric (JTJ. The term 
C/a 4 then behaves exactly as a radiation component. The main constraint 
on the amount of such extra radiation comes from nucleosynthesis, since 
the primordial abundances of light elements depend crucially on the balance 
between the expansion rate of the universe and the rate of microphysical 
reactions. The good agreement of observations with the predictions of stan- 
dard nucleosynthesis implies an upper bound on the number of nonstandard 
relativistic degrees of freedom in the universe and thus, in the context of 
brane cosmology, an upper bound on the value of the Weyl parameter C. 

Now, if the bulk is not strictly empty but contains a non-vanishing bulk 
energy-momentum tensor, then the Weyl parameter C is no longer necessarily 
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constant. In other words, the generalization of Birkhoff's theorem to brane 
cosmology no longer applies. This has been illustrated in particular in models 
with a bulk scalar field E] ■ 

In fact, in a realistic brane cosmological scenario, there is an unavoidable 
bulk component, which is simply due to the the gravitational waves, or bulk 
gravitons, produced by the brane matter fluctuations. Such bulk gravitons 
are created by the scattering of ordinary matter particles confined on the 
brane. The expected consequence of this phenomenon is to 'feed' the Weyl 
parameter by means of a transfer of energy from the brane into the bulk. 

This problem has been recently investigated in |7J |Hl El- Whereas the 
authors of [HI analyze generic forms of transfer of energy between brane and 
bulk using a four dimensional effective description, the case of pure graviton 
emission has been studied in HM [7j and LSR [8 . 

In HM, the authors consider a pure AdS bulk and analyze the generation 
of dark radiation by considering two distinct phases. In the low energy regime 
(p <C <y) the amount of generated dark radiation is obtained by directly 
equating it to the loss of energy density on the brane. In the high energy 
regime, the emitted gravitons are considered to remain gravitationally bound 
to the brane for the whole duration of this era, bouncing several times off it. 
To estimate the corresponding amount of dark radiation, the lost energy is 
reduced by a factor due to the energy dissipation resulting from the repetitive 
collisions with the brane. Because of the somehow crude distinction between 
the high and low energy eras, this analysis cannot follow the details of the 
transition between these two regimes. Moreover the authors of HM give their 
final result up to an uncertainty factor, estimated to be between 0.5 and 1, 
due to the multiple bouncing of gravitons off the brane. 

In LSR, a five-dimensional exact solution for the bulk is used, taking 
into account an energy flux from the brane into the bulk. The simplest 
solution of this type is the five-dimensional generalization of Vaidya's metric. 
In this way the brane trajectory, and thus its cosmological evolution, are 
determined self-consistently via the junction conditions, taking into account 
possible backreaction effects. The price to be paid for this is the necessity 
to make the assumption that the trajectories followed by the bulk gravitons 
are exactly perpendicular to the brane, assumption that is realistic only in 
the low energy regime. 

Although these two approaches are very different, they give estimates 
that agree within one order of magnitude. Moreover, these estimates are 
very close to the current observational bound on the number of extra rel- 
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ativistic degrees of freedom. Therefore, it is important to get an accurate 
determination of the amount of dark radiation expected for this model. In- 
deed, the improvement of observational constraints on the number of non 
standard relativistic degrees of freedom will allow us to put constraints on 
the cosmological evolution of the Randall-Sundrum model (possibly exclud- 
ing a long era of nonstandard H 2 oc p 2 cosmology) and on the parameter 
space of the model. A detailed analysis can also allow us to quantify the 
effects of the bounces of gravitons on the brane and to estimate to which 
extent backreaction effects can be relevant. The present analysis, finally, can 
be generalized to the case of models with extra matter in the bulk [5J G2] , 
eventually leading to stronger constraints. 

Our plan is the following. In the next section, we discuss the problem 
from an effective four-dimensional point of view. In section 3, we rederive 
the equations governing the cosmology of the brane and the trajectory of the 
brane in the bulk. Section 4 is devoted to the emission of bulk gravitons by 
brane particles. In section 5, we recall our model [H] based on the Vaidya 
metric. In section 6, we discuss the propagation of the gravitons in the 
bulk. In section 7, we present and discuss our numerical computations. We 
conclude in the final section. 

2 Effective approach 

As a starting point, we will discuss the cosmology of the brane and the 
influence of the bulk gravitons from a four- dimensional effective description. 
In order to do so, we will follow the approach of [TTj . 
We start with the 5-dimensional Einstein equations 

Rab ~ 7^9 abR + & B g AB = k 2 [T ab + S AB 5{y)} , (4) 

where the matter component, on the right hand side, consists of a bulk 
energy-momentum tensor Tab and distributional a brane energy-momentum 
tensor. The brane is located at y — 0, y corresponding to the proper coor- 
dinate normal to the brane. We assume that the extra dimension has a 
orbifold symmetry and that y = is a fixed point under this symmetry. 

In the cosmological extension of the Randall-Sundrum model |4j, Sab is 
the sum of a brane tension er, defined by 

K 2 a = V^6A = 6/1, (5) 
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and of the contribution of ordinary matter fields confined to the brane tab > 
i.e. 

Sab = —&hAB + tab , (6) 

where Kab is the metric induced on the brane. 

From the above five-dimensional Einstein's equation, using the Z2 sym- 
metry, one can derive ^T] effective four-dimensional equations that read 



(4) G 



with k\ = K 2 jj, and 
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where n" 4 is the unit vector pointing outward and normal to the brane. In 
addition to the usual four- dimensional matter energy- momentum tensor r^, 
three new terms appear on the right hand side of the effective Einstein equa- 
tions: the first one is quadratic in r Mi ,; the second one is the projection on 
the brane of the 5-dimensional Weyl tensor ^Cabcd] and the last one is the 
projected effect of the bulk energy-momentum tensor. 

Since we are interested here in homogeneous brane cosmology, all the ef- 
fective energy-momenta defined above have necessarily a perfect fluid struc- 
ture, i.e. 

rW = (p w +p w )«^+p (, V. ( 9 ) 
where is the timelike unit vector associated with comoving observers on 
the brane. 

It is then not difficult to show that Eq. (JJJ) gives a Friedmann equation 
on the brane of the form 



H 2 



P_ 
2a 



P + P 
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P 
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(10) 



where H = a /a is the Hubble parameter on the brane (a is the scale factor on 
the brane and the overdot stands for a derivative with respect to the cosmic 
proper time t). We have used the relation 
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Let us now discuss the various (non) conservation laws satisfied by the 
different energy density components defined above. We follow here the recent 
analysis of [10J, where the bulk energy- momentum tensor was associated to a 
five-dimensional scalar field. One can first show that the usual conservation 
equation for cosmological matter is modified into 

p + 3H (p + p) = 2T RS n R u s . (11) 

The right-hand side (evaluated at the brane position) represents the energy 
flux from the bulk into the brane. When the brane loses energy, as will be 
the case here via emission of gravitons, the right-hand side is negative. The 
factor 2 is a consequence of the Z2 symmetry. 

Moreover, the 4-dimensional Bianchi identities imply that the total energy 
density defined by p( tot ) = p + p^ + p^ + p^ satisfies 

p(tot) + 4Hp {tot) + #r (tot) % = 0. (12) 

Using = {p/<y){r% + 2p), = 2(k 2 / Kl)T AB n A n B , = and 

the (non) conservation equation (jllj) . it follows that the energy density for 
the dark component, by which we mean the sum of the components depending 
explicitly on the bulk, i.e. pr> = p^ + p^ w \ satisfies 

p D + 4Hp D = -2 (l + T AB u A n B - 2— T AB n A n B , (13) 
V a J p: 

where the right hand side is evaluated at the brane position. On the right- 
hand side of the above equation, one recognizes in the first term the energy 
flux from the brane into the bulk, —2 T A BU A n B . This means, not surprisingly, 
that the loss of energy inside the brane, will contribute to an increase of 
the amount of dark radiation. In the second term, the quantity T AB n A n B 
can be interpreted as the pressure transverse to the brane, due to the bulk 
component. In the case of a gas of gravitons, which we consider here, this 
pressure is positive and therefore this term tends to decrease the amount of 
dark radiation. The two terms on the right hand side have thus opposite 
effects. 
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3 Cosmology of a brane in AdS 

We now consider explicitly the bulk and study the trajectory of a brane with 
relativistic matter in a strictly AdS spacetime, with the metric 

ds 2 = -f(r)dT 2 + ^— + r 2 dx 2 , f(r) = ^ 2 r 2 . (14) 
f( r ) 

The trajectory of the brane can be represented in terms of its coordinates 
T(t) and r(t) as functions of the proper time t, which is also the cosmic time 
in the brane. The normalization of the velocity vector u A = (T, r, 0) then 
implies 

u = I , r, I . (15) 

One then needs the junction conditions for the brane, 

h c A S/ c n B = k 2 (t ab - \rh AB ^j , (16) 

where n A is the unit vector normal to the brane (pointing outwards) and is 
given by 

n A = ~ (7, VTT^, 0^ . (17) 



The spatial components of the junction equations yield the brane Friedmann 
equation, which can be expressed as 

H 2 p p 2 

— = 2 - J r-2- 18 

It is convenient to introduce the dimensionless quantities 

fjL a 

Note, using (jSJ) and that the brane tension can be expressed as 

a = 6fi 2 m 2 P , (20) 

where mp = 1/ y/El is the reduced Planck mass. We will drop the tildes from 
now on. 
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In terms of H, the components of the brane velocity and of the normal 
vector read respectively 

u A = + H 2 ,rH,0^j , n A = - (—^VT+hA. (21) 

The trajectory of the brane is obtained by solving 

dr f 2 H 2 V 2 P + P 2 / 00 N 
—— = — = r -. = r , 22 

dT T VI + H 2 1 + P 

where the second equality follows from ()21j) and the third from ()18j) . In the 
case of radiation domination, p = pi r -4 , where we have defined the radial 
coordinate r such that r = 1 for p = pi, where pi is the energy density at 
some fiducial initial time £«. To get the brane trajectory, one must integrate 



dr _ , pP^ + P 2 
dT~ T r* + p t ■ 

Explicit integration is possible in the high energy regime, pi/r 4 ^> 1, where 
one finds 

' -T + const. (24) 



r 



and in the low energy regime, where 



2 pi T + const. (25) 

In the regime of transition between the high energy phase and the low energy 
phase, i.e. for r ~ p { ^ 4 , one must resort to numerical integration. 



4 Emission of bulk gravitons 

We now consider the production of bulk gravitons by the cosmologically 
evolving brane. In particular, we wish to compute explicitly the components 
of the bulk energy-momentum tensor due to the gas of gravitons emitted by 
the brane. 

Let us first recall that the energy-momentum tensor due to a gas of mass- 
less particles is given by 

T AB = Jd 5 p6 ( PmP M ) V=9 IpaPb, (26) 
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where / = f(x A ,PA) is the distribution function of the gravitons (which is a 
scalar that depends on the spacetime position and on the momentum). 

Since the gravitons living in the bulk will be assumed to have been pro- 
duced only by emission from the brane, it will be convenient to start the 
computation of the energy-momentum tensor components on the trajectory 
of the brane. At the location of the brane, we can see the bulk gravitons 
from two perspectives. 

First, from the brane point of view, the bulk gravitons are seen as massive 
four- dimensional particles with mass m, three-momentum p and energy E, 
satisfying E = y/p 2 + m 2 . They are created by the scattering of two ordinary 
particles confined on the brane. The leading contribution to this process is 
given by the scattering if> if) — > graviton, where ip is a standard model particle. 
At the cosmological level, the production of gravitons results into an energy 
loss for ordinary matter, which can be expressed as 

T t + ^ + P) -lSf Cm ' (27) 

with 

i J% ^!^ |2/ ; (27r) 4 5(4)( } 

Ui 2 J (2tt) 3 2E 1 (2tt) 3 2E 2 ^ 1 1 JJK ! ^ F ^ 

(28) 

where M. is the scattering amplitude for the process in consideration, and 
the indices 1 and 2 correspond to the scattering particles (ip and i/j). 

Second, from the bulk point of view, the gravitons are massless particles, 
each graviton being characterized by a five- dimensional momentum, which 
can be decomposed into a spatial four-momentum and an energy, defined 
with respect to a reference frame. 

To make the connection between these two points of view, it is convenient 
to choose a reference frame associated to a (comoving) brane observer. Intro- 
ducing an ortho normal frame defined by u A , n A and e A , the five-dimensional 
momentum of any graviton can be decomposed into 

p A = Eu A + mn A + Pi4- (29) 

i 

The components along the vectors e A are denoted with a tilde to distinguish 
them from the components p l defined in the decomposition along the coor- 
dinate vectors (d/dx l ) A . If one substitutes this decomposition into (J26|) . one 



9 



gets for the mixed component 

T nu = u A n B T AB = J d 5 p5 (p M p M ) \f ZI g f PaPb n A u B = - J dm d 3 p y / . 

(30) 

By identifying the right hand sides of (fTTj) and (J2ZJ), one immediately finds, 
upon comparing ()28|) with the above expression, that the distribution func- 
tion at the brane location, for gravitons which are being emitted, is given 
by 

hem) K P) = ^ / |^ ^ E \M | 2 h h ^ (ft + P 2 - p) ■ 

2 m (2 n) J 2pi 2p 2 

(31) 

The summed squared amplitude is given by the expression 

Y.\M\ 2 = A^s> (32) 

where s is the Mandelstam invariant and 

A= -g s +g f + 4:g v , (33) 

where g s , gf and are respectively the scalar, fermion and vector relativistic 
degrees of freedom in thermal equilibrium on the brane. For the standard 
model, if all degrees of freedom are relativistic, one has g s = 4, gf = 90 and 
g v = 24. For simplicity, we will ignore the Bose or Fermi corrections in the 
distribution functions and assume that the particles on the brane are char- 
acterized by a Boltzmann distribution of temperature T. The distribution 
function at emission, f( em ), can then be computed explicitly and one finds 
(neglecting the masses of the scattering particles) 

ke m) (m, P) = ^E- 2 - 3 e-V^/T . (34) 

We now have all the elements to compute all the components of the bulk 
energy-momentum tensor due to the emitted gravitons at the brane location. 
They are given by 

Tjr> = /,i m ^p| /(e ,„ ) = JL^r» = ^^V, (35) 
r<r» = -/*^ 3 pf/, em , = -f^ K 2 r« = -^^V, (36) 



rir> = J dmd ^ hm) = l±^T^ e §± K y, (37) 
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with E = VP 2 + iTT'' 2 as stated before. In the second equalities, we have 
replaced the temperature by the energy density, using p = (7r 2 /30)g*T 4 , 
where g* is the effective number of relativistic degrees of freedom. The above 
components can be interpreted physically respectively as the energy density, 
the energy flux and the lateral pressure of the emitted bulk gravitons as 
measured by an observer at rest with respect to the brane. 

5 Cosmological evolution in the Vaidya model 

In this section, we deviate from the preceding analysis and consider the 
evolution of the dark radiation in the context of the Vaidya model introduced 
in [Hj (see also The interest of the Vaidya model is to work with an 

explicit five-dimensional realization of the effective equations introduced in 
section 2, where there is a non zero energy transfer between the brane and 
the bulk. This however imposes that all the gravitons have to be assumed 
to be emitted radially in the five-dimensional bulk so that the bulk energy- 
momentum tensor is of the form 

Tab = cr B k A k B , (38) 

where k A is a null vector. Assuming spherical symmetry, Einstein's equations 
(@J) with such an energy-momentum can be solved analytically: this is the 
Vaidya solution jT3|, ordinarily used to describe a radiating relativistic star. 
Its metric is given by 

ds 2 = -f( r ,v) dv 2 + 2drdv + r 2 dx 2 , (39) 

with 

/M) = A 2 -^- (40) 

For a constant C, one recovers, after changing the light-like coordinate v into 
the static time coordinate t, the familiar five-dimensional AdS-Sch metric. 

It is instructive to apply the effective equations derived earlier in this 
particular context. First, one can normalize the null vector k A such that 
kj\u A = 1, in which case the projections of the bulk energy- momentum tensor 
f)38p are given by T nu = —a B and T nn = a B . As a consequence the non- 
conservation equation (jTT| for the energy density on the brane reads 

p + 3H {p + p) = -2a B , (41) 
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and Ob, up to the factor of 2 due to the Z 2 symmetry, directly represents 
the energy loss in the brane due to the production of bulk gravitons. The 
dependence of Ob is given explicitly in (|36jl. and in terms of the brane energy 
density, 

ob oc p 2 . (42) 

The Weyl parameter C can be related to the "dark" component energy 
density defined earlier via 

Pd = (43) 

Substituting in the evolution equation (JT3J), one gets (still using the implicit 
tilded quantities defined in ()19|l) 



C 

a 



, - -2 (1 + p) T AB u A n B - 2HT AB n A n B . (44) 



With the explicit bulk energy-momentum tensor (jSHJ)? T nu = —<tb and T nn = 
Ob, Eq. (JH|) reduces to 

^ = 2a B (l + p-H) , (45) 

where one recognizes the result already derived in [S], but there from Ein- 
stein's equations rather than from the effective equations like here. The 
Friedmann equation, when one neglects the dark radiation contribution, gives 
during the high energy regime 

H~p + 1-— . (46) 
4p 

This shows that there is a remarkably precise compensation, at leading order 
and next to leading order, between the energy flux and the transverse pressure 
in the dark radiation equation. And the production of C is governed by 



a 4 ~ 2p 



i - . « p. (47) 



As mentioned above, the drawback of the Vaidya description is the as- 
sumption that all gravitons are radial. As seen in the previous section, the 
distribution of emitted gravitons is not radial and T nn + T un is not zero. 
Therefore, if one substitutes the explicit values of T nn and T un in (jUj) we no 
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longer get the very precise compensation observed in the Vaidya case and the 
global sign is positive (since 315/(2 10 7r 3 ) > 3/(47r 4 )), so that the production 
of dark radiation seems to be driven in the high energy regime by a term on 
the right hand side proportional to p 3 , rather than proportional to p as in 
the Vaidya description. Obviously, this would result in an enormous amount 
of dark radiation, far above the estimate of jS] , and this would ruin the sim- 
plest brane cosmology scenario, since the estimate of jS] was barely within 
the nucleosynthesis bounds. 

The above analysis is however incomplete. In the Vaidya description, the 
gravitons are emitted radially inwards and are therefore lost for the brane 
once there are emitted. Some of the non radial gravitons, however, can come 
back onto the brane after their emission, as it will be shown in the next 
section, and thus influence once more the evolution of the Weyl parameter. 
Because of the Z2 symmetry, these gravitons will be reflected by the brane 
(we ignore here the decay of a bulk graviton into brane particles) and will 
contribute only to the transverse pressure term. The effect of these old 
gravitons will thus be to reduce the amount of dark radiation that would be 
computed naively by considering only gravitons being emitted. 

6 Graviton trajectories in the bulk 

After the gravitons are emitted, they move freely in the bulk, each individual 
graviton following a null geodesic. The null geodesies in five-dimensional 
AdS were studied in |14j . We summarize here the main results. Using the 
symmetries of the metric (|14jh one can identify the first integrals for the 
geodesic motion 



where A is any affine parameter. For any graviton one can choose the affine 
parameter so that the tangent vector of the null trajectory is identified with 
the physical momentum p A . Introducing the notation p T = rp T and p r = 
p r /r, the above first integrals become 



f(r) 



dT 
~d\ 



£ 




(48) 



£ = rp 



(49) 



and 



V = rp. 



(50) 
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Using these conservation laws, it is easy to determine the trajectory of 
the gravitons in the bulk spacetime. In order to do so, one can compute 



dr _(dr\ I dT 
~dT ~ \d\) 1 \d\) ' 

Since p A pA = 0, one gets 



/ hr ■ ( 51 ) 



(pry = s 2 - V 2 , (52) 

which means that p r is also a constant of motion. It is then useful to introduce 
the parameter 

V = f , (53) 

which will be constant along the null geodesic trajectory, and in terms of 
which the trajectory in the time/extra-dimension submanifold is given by 

*=W. (54) 
Integration gives the graviton trajectory 

V \r r*/ 

It is also useful to relate the bulk-based description used above with the 
brane-based approach introduced earlier. Using the decomposition (}2T?j) with 
the velocity and normal vectors given in (}2T]) . one finds 



f = y/i + H 2 E - Hm, f = HE - VI + H 2 m. (56) 
It is immediate to invert this system, to get 



m 



Hf- Vl + H 2 p r , E = Vl + H 2 p T - Hp r . (57) 



For some values of the parameters, the trajectory of the graviton, leaving 
the brane at the emission point, can cross again the brane trajectory. It will 
then be reflected by the brane, as a consequence of the Z2 symmetry. 

It is not difficult to compute the new five-dimensional momentum of the 
graviton after the reflection on the brane, by considering the reflection from 
the brane point of view. The momentum parallel to the brane p is conserved, 
as well as the energy E. Only the transverse momentum is affected and sim- 
ply changes its sign. The momentum along the spatial direction orthogonal 
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to the brane is embodied by the mass m. We will adopt the convention that 
m is positive when the momentum is outwards, with respect to the brane, 
and negative otherwise. The reflection of the graviton by the brane is thus 
governed by the simple laws, 

E — > E, p — > p, m — > —m. (58) 

In the bulk-based point of view, this translates into 

f out = (l + 2H 2 )g n -2HVTTlPpr n , (59) 
p r out = 2HVT+lPpf n -(l + 2H 2 )pl n . (60) 

We can use the previous results to infer the evolution of the graviton 
distribution function throughout the bulk. In principle, it is governed by the 
Liouville equation, which, in general relativity, reads 

P ^ + rS BP V^ = o. (6i) 

However, this equation is in fact no more than the implementation at the 
level of the distribution function of the fact that each individual particle 
follows a geodesic. To solve explicitly the Liouville equation it is therefore 
simpler to use directly the solutions for the geodesic trajectories. 

This enables us to write the emission distribution function, at the brane 
location, in terms of the bulk-based momenta, by simply substituting the 
above expressions in (|34j). Moreover, using the constants of motion along the 
graviton geodesies established above, one can deduce the expression for the 
distribution function off the brane as well. One gets 

ry ry ry 

r,p T ,f, p) = /(em) (Tem, a em , f , f , p), (62) 

"em "em "em 

where the time and brane scale factor at emission depend on t, r, p T and p r 
and are obtained by tracing the graviton geodesic back in time until it crosses 
the brane trajectory. The rescaling of the momenta on the right hand side are 
just a consequence of the conservation laws (f4"H|) and (jHOj) (the same rescaling 
for p r follows from the normalization of the five-dimensional momentum). 
By construction, the above distribution function must be a solution of the 
Liouville equation (|6*T|) in the bulk, as one can explicitly check. 

The above expression describes gravitons which have been emitted by the 
brane and have not returned back onto the brane. But, as already mentioned, 
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the situation is in fact more complicated because some of the gravitons that 
have been emitted by the brane can come back onto the brane and be reflected 
back into the bulk with a different momentum. Therefore, at each point along 
the brane trajectory, one must add to the gravitons that are being emitted all 
the older gravitons that are being reflected by the brane at the same instant. 

Actually, the majority of gravitons emitted in the high energy regime will 
be reflected at least once by the brane. In the period in which H 2 oc p 2 
(i.e. H 3> 1) the brane is indeed moving relativistically with respect to the 
frame defined by the metric (14). The gravitons that are not emitted exactly 
orthogonal to the brane will see their radial momentum boosted, and will 
thus move in this frame in the same direction as the brane. This can be seen 
by inspection of eq. (56): for H ^> 1 and E not too close to m, p r is positive. 
As H decreases, the brane will eventually slow down, and the gravitons will 
bounce off it. 

Let us consider a fiducial time, characterized by to, in the history of the 
brane. In order to compute the increment of the Weyl parameter, via (|44|) . 
one must add to the newborn gravitons, which are just being emitted, the old 
gravitons that are being reflected by the brane. For these old gravitons, not 
surprisingly, the net contribution to the energy flux T nu vanishes, since they 
are just reflected, neither absorbed nor created, whereas the incoming and 
outgoing contributions to the pressure are equal in magnitude and opposite in 
sign. For these gravitons, therefore, we need to compute only the component 

?l n) = Jdm£p^f im) . (63) 

Let us now write explicitly the contribution to /j„ from gravitons which 
make at the time to their first bounce since their emission. It will be conve- 
nient to introduce the following parametrization 

H = sinh (3, (64) 

for the Hubble rate, and for each graviton, the parameter x such that 

E = pcosh.x, m = psmhx, (65) 

where p = |p| is the norm of the three- momentum parallel to the brane. 
Using the results obtained previously, one finds that the energy E% and mass 
mi measured in the brane frame at the time of emission t\ are given by the 
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Figure 1: Trajectory of the brane and typical trajectories of gravitons in the 
metric (14). The brane trajectory is the solid thick line, the steepest part 
on the left corresponding to the high energy regime p> 1. The dashed line 
describes a graviton produced at t 2 , bouncing off the brane at t\ and again 
at to- The dashed-dotted line corresponds to a second graviton, that, after 
being reflected once by the brane, falls into the bulk. 



simple expressions 

Ei = —p cosh (Pi — P + x) , mi = —p sinh (Pi — p + x) , (66) 
a\ ai 

where 

Pi = Pi(Po,x), (67) 

is determined by tracing backwards the null geodesic (see Fig. 1) that in- 
tersects the brane trajectory at t until the previous intersection. More ex- 
plicitly, Pi is obtained by solving the equation T b (r) = T g (r), where T b (r) 
gives the time T as a function of the position r of the brane through eq. (|23|) . 
whereas T g (r) is given by the graviton geodesic equation (|H5|) . 
Substituting 

{1) hn) = ^E^mle-^ T \ (68) 
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in the integral (|63|). one finds, after integration over p, 



H) r iin) ( a \ 71 A K 2 f sinh 2 a;sinh 3 (A(x) - O + x) / a \ 5 8 

In a similar way, one can compute the fraction ^/m due to the gravitons 
for which this is the second bounce on the brane since their emission. Con- 
sider such a graviton, which was emitted at time t 2 , was reflected by the brane 
once at time t\ and comes back again onto the brane at time to (t 2 < t\ < t ). 
Just before the first bounce, the extra-dimensional momentum is given by 

mi" 1 ' = pi sinhxi = — ni]™^ = — — psinh (/3i — (3 + x) . (70) 

Hence, at emission, we had 

m (out) = ^o pginh (/?2 _ 2/3i + (3o _ x y (7i) 
a 2 

and thus 

E {<rut) = ao p cogh ^ _ ^ +/3o _ x y (72) 
a 2 

Consequently, the contribution to the transverse pressure is of the form 

(2) T (in) (R \_ 7]A 2 f > sinh 2 x sinh 3 (/j 2 (x) - 2/?i(s) + (3 - x) 
nn [f 0) " 29 vr 4 K J cosh 8 (/? 2 (x) - 2/3x(x) +/3q-x) 

T|(x), (73) 



\ -5 

a 



\a 2 (x) / 

where we have stressed the dependence on x in the integrand (there is also 
a dependence on /3o). It is then straightforward to generalize to the contri- 
bution of gravitons with any number of bounces. 



7 Numerical results 
7.1 Source terms 

Whereas the contributions from the newborn gravitons to the energy flux 
and pressure can be computed analytically, one must resort to numerical 
integration to compute the contribution from old gravitons. For the latter 
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case, we need only the contribution to the pressure, since the contribution 
to the energy flux vanishes as explained earlier. Therefore, one can write the 
evolution equation for the dark component as 



p D + AHp D 



-p{em) _ jy(b) 



(74) 



with 



V (em) = _ 2 (i + p ) ri em ) - 2HT^ m) , V ih) = 4H1$p . 



(75) 



The above equation can be rewritten as 



d 

da 



(a* PD ) = a -V^ - a -V^ 



(76) 



where p( em ) can be computed analytically by means of f)36H37j) . 

We have computed numerically the function V^ b \a) using the general- 
ization of Eq. (J75|) . taking into account gravitons that have made up to 10 3 
bounces before hitting the brane when its scale factor was a. We have checked 
that neglecting the effect of gravitons that have made more than 10 3 bounces 
does not change appreciably our results. We have also assumed that the bulk 
contains no gravitons at the initial time ti. 

Numerically reliable results could be obtained for values up to pi ~ 10 3 . 
Taking the lowest value of \x compatible with the data (from small-scale 
gravity experiments), pT x ~ 0.1 mm, one gets that the brane tension has 
to be at least of the order of the TeV, which corresponds to a fundamental 
Planck scale M 5 = «r 2 / 3 ~ 10 8 GeV. As a consequence, the highest initial 
energy density which makes sense is pi ~ M^/a ~ 10 20 . 

In order to compare, at each instant in the history of the brane, the 
contribution directly due to the emission of gravitons with the one due to 
the reflection of older gravitons, we have plotted in Fig. 2 the two terms on 
the right hand side of eq. (|75jl . One observes that at early times the dominant 
contribution is that due to the emission of gravitons, which is not surprising 
since the bulk is assumed to be empty initially. The bulk is then gradually 
filled with gravitons and some of them can come back onto the brane and 
be reflected. They contribute to the transverse pressure effect, which can 
be seen on the plot to build very quickly. In the intermediate phase, the 
pressure effect dominates the emission effect so that the net source term 
for the dark radiation is negative. However, the two integrated effects are 
very close in amplitude, which means that the compensation observed in the 
simple Vaydia model is still working in this case. This is also the origin of 
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Figure 2: Contributions to the emission of dark radiation (for pi = 1000) 
from graviton emission and from the pressure of bulk gravitons. V^ b > and 
are defined in eq. (J75|) . 



the numerical difficulties of the present analysis: the result we are looking for 
is the small difference of two very large numbers. One can indeed see that 
both terms on the right hand side of eq. (J75J) scale roughly as pj, whereas 
their difference scales approximately as p^. Notice that, for this reason, the 
quantities on the vertical axes of figures 2 and 3 have been rescaled by a 
factor p\. 

It is also important to stress the necessity to take into account the multiple 
reflections of gravitons on the brane to get a correct evaluation of the total 
effect. To illustrate this point, we have plotted in Fig. 3 the contribution 
from gravitons for which this is the first bounce in comparison with the 
cumulative contribution from gravitons that have made up to 10 3 bounces. 

7.2 Dark radiation 

Our main goal is to compute the dark radiation globally produced in the 
process. At very low energy, dark radiation is produced at a negligible rate, 
so that one can consider that there is an asymptotic constant value for the 
Weyl parameter C. In the end, this asymptotic value for C depends only 
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Figure 3: Contribution to the source term (a 3 /H)V^ from gravitons that 
have been reflected only once by the brane, with respect to the total term 
(numerically, up to 10 3 bounces). 




on the initial energy density in the brane and on the number of relativistic 
degrees of freedom during the high energy phase. It will be convenient to 
express the final amount of dark radiation as its ratio with respect to standard 
radiation energy density 

e D = ^ (77) 

Prad 

Such quantity is constrained by cosmological observations. The amount 
of nonstandard radiation in the early universe is usually measured in units 
of extra neutrino species AN V . The conversion factor between AN V and 

1 /3 

e D is AN U = (43/7) (gT cl /g*) e D , where gT d = 10.75 is the number of 
relativistic degrees of freedom just before nucleosynthesis, more precisely just 
before the electron-positon annihilation. The factor (g" uc / 'g*) 1 ^ accounts for 
the change in the number of relativistic degrees of freedom since the era in 
which dark radiation has been produced. Assuming g* = 106.75, this gives 
en — 0.35 AN U . Constraints on AN U depend upon the kind of observation 
one takes into account, but a typical order of magnitude is AN V < 0.2 fE^ lTo]. 
which gives an upper limit < 0.07. According to both the analyses [151 ITd]. 
the favored value of AN U turns out to be negative. Although this possibility is 
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allowed by a negative value of the Weyl parameter C, this would correspond 
for the metric (JTJ to a naked singularity in the bulk, a situation that is 
certainly not appealing from the theoretical point of view. 

In the evaluation of the total amount of dark radiation produced, an 
important quantity is the energy loss due to production of gravitons, which 
depends on the number of relativistic degrees of freedom confined on the 
brane. In LSR, the expression for the energy loss is given by 

ob = ^ 2 p 2 , (78) 

where 



212625 , . s , , . g 

with 



« = ^n^C(9/2)C(7/2)^ (79) 



g = ({2/3)g s + 4g v + (l - 2~ 7 / 2 ) (l - 2^l 2 ) g f ) , g* = g, + g v + (7/8)g f 

(80) 

This is in exact agreement with the revised version J7| of HM. 

In the present work , o B = —Tj?™) and the expression for 7^ m ) is given 
in so that a, defined as in (|78|). is given by 

212625 A 
647T ' g% 

The difference between the two above values for a comes from the fact that 
Fermi-Dirac or Bose-Einstein distribution functions were used in HM and 
LSR, whereas here we have treated all the particles on the same footing 
and assume simply a Boltzmann distribution. To make a precise comparison 
between the previous estimates and our numerical results, we will use the 
analytical estimates of HM and LSR with the value for a given in (|8ip. For 
the particle content of the standard model, the values of a as given by eq. (|79j) 
and (|8ip differ only by about the 5%. 

Let us now recall these previous estimates. In HM, the amount of dark 
radiation produced in the high energy regime can be expressed as 

e D ~F-\n{ Pi /2) : (82) 

where T is an efficiency factor (denoted a in HM) with 57r/32 < T < 1. 
The amount of dark radiation produced during the whole low energy regime 
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(assumed in HM to start at p = 2) is a/2. Therefore in our comparison 
we will take for the total fraction of dark radiation as estimated in HM the 
expression 

e* M = ^[2 + Fhi(p i /2)] , Pi » 2 . (83) 

The value of en as given by LSR is obtained by solving a system of 
differential equations including ()45|) . In the limit of large pi one gets eo —* 
a/A. 

Let us now evaluate the amount of dark radiation produced in our nu- 
merical approach. Starting from the equation 

p D + AHp D = V {em) - V {b) , (84) 

where the source term is given explicitly above, one finds 

£fl(t0 ) = fS^L = I f" dt ( a M]' « _ pW) . (85 ) 

Prad{t ) Pi Jti \ Oi J V ' 

We have integrated numerically this equation and plotted in Fig. 4 the 
evolution of the ratio €d = Pd/p&s& function of the scale factor. Remarkably, 
for large enough initial values of the energy density on the brane pi, the dark 
radiation component can dominate the brane matter energy density at early 
times. This however does not seem to invalidate the implicit assumption that 
one can neglect the backreaction of bulk gravitons on the bulk geometry. 
Indeed, the dark energy domination, pr> (a) > p(a), is effective only in the 
high energy regime p ^> 1. And in this regime, the cosmological expansion 
is dominated by p 2 , which remains much larger than pp. Nevertheless, it 
might be interesting to explore the consequences of the fact that in the very 
early evolution of the Randall-Sundrum universe, most of the energy density 
of the universe should be in the form of dark radiation (or bulk gravitons) 
gravitationally bound to the brane. 

In Fig. 5, we compare the numerical estimates of the present work with 
the analytical estimates of HM and LSR. We took for a the expression (|81|) 
with all the degrees of freedom of the standard model. The curve correspond- 
ing to HM comes from eq. ()83|) with the lowest value of T . 

Due to the difficulties mentioned above, the range of our numerical anal- 
ysis is limited to values of pi smaller than about 2 x 10 3 . In this range, we 
get an accurate description of the behavior of as a function of p^. We 
see that is a slowly increasing function of p.;. At low values of pi, the 
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Figure 4: Evolution of Ed = Pd/p for different values of pi. 



numerical results are close to the estimate of LSR . Indeed, for p { < 1, the 
effect of graviton bounces can be neglected, and the Vaidya description an- 
alyzed in section 5 gives a good approximation, whose results agree, in the 
limit pi <C 1, also with those of HM. Going to higher values of p$, in the 
regime of validity of our analysis the value of en is below the lowest bound 
estimated by HM, but it gets closer and closer as pi increases and one can 
expect it to become higher already for pi of the order of few thousands (re- 
member that the maximal value of pi compatible with constraints on the 
Randall-Sundrum model is pi ~ 10 20 ). It is remarkable that already for the 
(relatively) low value pi ~ 1800, tr> > 0.02, which is not far from the upper 
bound imposed by CMB and BBN observations [TBI llb'|. 



8 Conclusions 

In the present work, we have computed the amount of dark radiation pro- 
duced during the cosmological evolution of our brane-universe. Previous 
estimates were based on relatively crude approximations: in one case, the 
evolution of the brane was separated into a high energy phase and a low en- 
ergy phase; in the other case, all bulk gravitons were supposed to be radial. 
However, both estimates were pointing to values very close to the current 
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Figure 5: Comparison of the numerical results of the present work with the 
estimates of HM (upper curve) and of LSR (lower curve) for the amount of 
dark radiation eo = Pd/p as a function of the initial energy density on the 
brane pi. 



bound on extra relativistic degrees of freedom, which was motivating a more 
detailed study of the question. 

To go beyond the previous approximations, the present work uses a nu- 
merical approach, which enables us to treat smoothly the transition between 
the high and low energy regimes and to deal with non radial gravitons. Our 
initial objective to compute precisely the amount of dark radiation produced 
is however hampered by a problem of numerical precision, and we have been 
able to do this computation only for moderate values of the initial brane 
energy density. The reason for this limitation lies in a remarkable compen- 
sation between two opposite effects: the emission of bulk gravitons by the 
brane, which contributes positively to the dark radiation, and the pressure 
of old gravitons bouncing off the brane, which contributes negatively to the 
dark radiation. Numerically, we have been obliged to compute these two 
effects separately, and the net effect, in which we are interested, comes from 
the difference of two quasi-equal huge numbers, which is difficult to control 
numerically. 

With the present analysis, we could estimate the evolution of the dark 
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radiation component as a function of time, finding that in the very early 
stages of a radiation dominated brane universe, energy as dark radiation 
could easily exceed the amount of energy in brane matter. 

In the range of initial densities where we can rely on our numerical compu- 
tation, we have been able to compare our result with the analytical estimates 
obtained previously and our results agree with their order of magnitude. 

The increase of the amount of dark radiation produced with the initial 
energy density on the brane seems to indicate that a significant amount of 
dark radiation will be produced for an extended high energy (non standard) 
period in the early universe. However, we cannot safely extrapolate the 
present analysis to higher values of the initial energy density on the brane 
and we believe that an extension of our numerical computation further into 
the very high energy regime will tell us whether a long period of nonstandard 
cosmology is going to be soon ruled out by observations. 
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